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Abstract 

Transitions to measure synchronization in two coupled 0^ lattices are inves- 
tigated based on numerical simulations. The relationship between measure 
synchronization (MS), phase locking and system's total energy is studied both 
for periodic and chaotic states. Two different scalings are discovered during 
the process to MS according to phase locking. Random walk like phase syn- 
chronization in chaotic measure synchronization is found, and phase locking 
interrupted by phase slips irregularly is also investigated. Meanwhile, related 
analysis is qualitative given to explain this phenomenon. 
PACS numbers: 05.45.Xt, 05.45.Pq, 05.45. Ac, 05.45.-a 
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The phenomenon of synchronization has been studied since 17th [0], but the related 
research never stopped. Recently, in order to explore the collective behavior in complex 
systems, this topic has been paied more attention than ever before and became a hot focus 
in many science fields 0. Although synchronization is first mentioned in periodic systems, 
synchronization in chaotic systems attracted more attentions due to its rich dynamical prop- 
erties and potential applications in practical systems 0. During the research process, the 
original concept has been largely extended, many new phenomena had been discovered, 
phase synchronization @], lag synchronization [||, generalize synchronization |^ had been 
the familiar words in new related works. Till to now, most attentions had been concen- 
trated on dissipative systems, the related research for conservative system, the hamiltonian 
system, has been neglected and only few works appeared 0. Hamiltonian system stands 
as the central topic of both classical and quantum mechanics, and a large variety of practi- 
cal systems can be well approximated by Hamiltonian formalism even at weak dissipation. 
Thus, it is important to understand the behavior of Hamiltonian systems in their possible 
synchronization process. 

The former works abut Hamiltonian systems are mainly about their collective behav- 
ior in globally coupled high dimension systems or studied their transitions to equipartition 
from non-equilibrium state based on statistical methods few paper had concerned to 
investigate their detail dynamical behaviors during the system's evolutions. The reason 
mainly because that Liouvill's theory request the phase space be conservative, it prevents 
the full collapse of the orbits which related with complete synchronization and phase syn- 
chronization, this means the classical concept for synchronization developed in dissipative 
system is unsuitable for this kind of system. But this situation has been changed in 1999, 
a new phenomenon, measure synchronization (MS), for Hamiltonian systems is discovered 
by Hampton and Zenette 0. In their work, MS is defined as all orbits cover the same 
region in the phase space with identical invariant measures, and the transition to this state 
is characterized by a critical coupling strength. Although MS has been found and critical 
coupling is defined, that are still many questions still leaving unsolved. For example, the 
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route from desynchronization to MS is unclear, and how about system's variables behave 
and what's the innate character for MS still need to be cleared. Also, the relation between 
chaos and MS need to be considered again, especially, the dynamical behavior for MS in 
chaotic state is also valuable to study carefully. In this paper, we will concerned mainly on 
these questions, and try to explore the way from non-MS to MS in coupled Hamiltonian 
systems. 

Let's consider two coupled classical (^^ lattices, the system's Hamiltonian reads 



= + ^ + (1) 

the direct physical image is two coupled point masses each of them oscillated in an quartic 
potential. Due to the quartic potential, system could displays complex behavior and may 
express new property out of intuition. By changing coupling strength e and initial values, 
both period and chaos can be realized. The corresponding canonical equations are 

^1 =Pi, Pi= -ql + e{q2-qi) 

^2 =P2, P2=-ql + e{qi-q2) (2) 
In this model, the system's total energy is conserved, E{t) = Pi^^^^+^aW _|_ ii W+g2 W _^£(^q^(^f'^ _ 
q2{t))'^ = E{0) =Const, and initial values act as an independent additional parameters. If we 
set each lattice's initial positions at zero, then total energy will be independend to coupling 
strength, this will be used in the following. 

Once two lattices reached MS, then any macroscopic values will equal if only the average 
time long enough. One example is the energy of each lattice, defined as 

e. = ^ Jjpm/^ + qtit)/A]dt, z = 1, 2 (3) 

before MS, two orbits occupy different domains which never intersect, of course ei ^ 62, 
this situation continued till MS reached, after the critical coupling, two lattice's energy be 
equaled due to their identical invariant measures. In Fig. 1 we choose initial condition 

qi = 0, q2 = 0, 

pi = 0.1, P2 = 0.2 (4) 



and plotted each lattice's energy Cj versus e. During our numerical simulations, third order 



bilateral symplectic algorithm [^] and small time step h = 0.001 are used to keep the 
system's conservative property. The critical coupling, Sc = 0.0032, which separate non-MS 
and MS is clearly shown in this figure. Meanwhile, we also plotted two lattice's total energy 
e = Ci + 62 versus e in the same figure for comparison. It is shown that energy e has its 
minimum on the critical point Ec, below and above this point e decreased and increased 
monotonically. This made us reminisced phase synchronization which has been warmly 
studied recently [|] |^ ||l3l, made it interesting to investigate the phase locking, or named 



phase synchronization, in this system. In Hampton's work, the phase difference = 61 — 62, 
with ^12 denotes each orbit's angle and defined on the real line rather than on the circle, 
between two orbits is found to obey different rules. Before MS, A6 increases linearly with 
time goes on, and above MS, the difference oscillates with a characteristic frequency which 
around a constant value. But only these can't explain the minimum, the deep reason need 
to be explored. 

More information about MS and phase locking can be found in Fig. 2, in Fig. 2(a) and 
(b) we adopted e = 0.002 < Ec, two orbits stay in different domain, and phase difference 
increases linearly with speed rj = 0.1. While e = 0.0035 > Ec, MS reached and phase 
difference is oscillates within {—it, it) with frequency / = 2.12, as what is shown in Fig. 2(c) 
and (d). In order to explore the relation between phase difference and coupling strength, we 
plot f], which is defined as the speed of phase difference between two orbits before MS, and 
/, which is defined as the oscillate frequency of phase difference after MS, versus e in Fig. 
2(e) and (f). It is observed that, before MS, rj decreases gradually as e increases and trend 
to zero on Ec, while / increases gradually from zero as e increases from Ec, this is coincide 
to Ref. §. 

During our numerical simulations, we found that phase locking exists two different trends 
corresponding to non-MS and MS respectively, future study shows that its right the reason 
why the energy e show its minimum at Ec- Before the critical point, MS doesn't exist and 
orbits have different main frequencies, but we found that: except 77 decreases as e increases. 



the distribution of phase difference </? (model with 27r) has a obviously trend. In Fig. 
3(a) we plotted the distributions for ip under three different coupling which all below the 
critical point. From this figure, we can find that as e increases, the probability for = tt 
changes to higher, we remark here that tt is the biggest phase difference for two orbits after 
modeled di and 82 with 27r. In Fig. 3(b) we also plotted the probabihty distributions for 
distance difference Ag = |gi — ^2!, it is found that with (p trends to tt in distribution, the 
probability for large Agr increases as well. After e > £c, as coupling increases, the oscillate 
frequency rise and the amplitude decrease, meanwhile, the distribution for Ag contracted 
to small scope. 

Already had these ideas in mind, the reason for the minimum of energy e on the critical 
point is straightforward. As the coupling increases, the relation between two phases become 
stronger, although they are unlocked, but they have the trend to reach phase locking, this 
can be found from 77 decreases as e increase. As we already known, for Hamiltonian system, 
the classical phase synchronization can't be reached due to its conservative character, but 
as coupling increases, two orbits trend to reach it. On the way to phase synchronization, 
the first A^ should be tt, because tt stands for the largest phase difference if we model phase 
by 27r. This directly made the probability of large Ag increased and thus induced energy e 
decrease in this process. Advanced increase the coupling, e > Sc, then two orbit become MS 
and be phase locked, under this situation, the coupling strength is strong enough to reach 
MS, but is weaker to bind the phase oscillate in a small amplitude, continue increase £, the 
amplitude will drops from large value tt, where e slightly larger than Ec, to zero. During this 
process, phases are locked and the oscillate amplitude decreased, induced the distribution 
for Aq inclined to short distance, at last made energy e tend upward. In order to give a 
more convincing expression, we plot the probability of tt (denote by p^) versus e in Fig. 
4(a), log(p^) for vertical axis and log(£c — s) for horizontal axis. A fine scaling is found 

with i/ = 1.96 ft! 2.0 and x denotes as constant. Also, in Fig. 4(b) we plotted amphtude A 



versus e in MS region, this time the scahng is found to be a exponential mode 



(6) 



with a = 1.0 X 10~^ and Aq, B constants. 

Under the situation mentioned above, Sc not only stands as the critical point for MS and 
the energy minimum, but also stands as the critical point for two different trends of phase 
locking, i.e., before MS, two orbits are trend to achieve the first possible locking = tt, 
after MS, the amplitude of oscillation dropped from tt and trend to approach zero. But 
we should also noticed here that all what we have discussed is just within period, no chaos 
appeared, so a natural question is how about these phenomena within chaos, does different 
phase locking trend still exist, and what's the new characters for MS when system is chaos. 
All these questions aroused our interests to investigate MS in coupled chaotic Hamiltonoian 
systems. 

In order to get chaos, this time we choose initial condition 



Similar to Fig. 1, we plotted energy Cj and total energy e versus coupling e in Fig. 5(a), 
different to the above result, energy behaves a disordered way in chaotic regions, made it 
difficult to get the critical point for MS. Through tracing the time evolution of their orbits, 
we found that the uneven behavior of is closely related to chaos, in other word, its just the 
conjoint action between chaos and MS that lead to this irregular behavior. On one hand, the 
definition of MS requests each lattice's energy be equaled, just as what happened in periodic 
state, but on the other hand, the internal stochastic character keep their orbits moving in 
an irregular way. Thus, ei — e2 can be reached only under limit T ^ oo. However, we 
found the behavior of summation e is more regular than e,, despite of some abrupt windows. 
Through calculate system's largest Lyapunov exponent A, it is found that these windows are 
corresponding to period which embed in chaos. In Fig. 5(b) we plotted e and A versus e in 



gi = 0, ^2 = 0, 



Pi = 0.1, p2 = 0.4 



(7) 
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the same figure, compare with Fig. 5(a), it is shown that e decreases more smoothly then 
Cj in chaotic regions. This result is somewhat similar to the behavior of e in region e < Sc 
plotted in Fig. 1. As we have known, the trends of phase locking plays important role for 
the minimum of e in periodic state, so we may ask does these the same to chaos, if not, 
what's the reason made e decrease in chaotic state and how about the phase locking now? 

On important difference between this work and Hampton's is, in our model, chaos is 
closely related with MS, and the largest Lyapunov exponent is a suitable order parameter 
for MS. We choose e = 0.025 > 0.012 = Ec, with initial condition Eqs. 7, system is chaos 
with A = 0.08. Phase difference A^^ versus time t under this coupling strength is plotted in 
Fig. 6. It is shown that A6 behaves a random walk process around zero, and there also exist 
some platforms with different wide. The inset figure is enlargement of the rectangle region 
in original figure, there the details about the platform is shown. In fact, each platform is 
in the state of phase locking, in these regions, although A6 oscillate all the time, but its 
amplitude is small than vr. We also found that even in regions where phase are unlocked, the 
phase difference shows a 27r like slips. The same phenomenon also can be found in others 
chaotic states with different e > Ec, the only difference is, as e increases, the random walk 
more near zero, and more time two orbits be phase locked. 

The 27T slips in phase locking made us reminisced the phase slips in phase synchronization 



of dissipative systems [|T2| [|T^ where various models were used and rich phenomena were 
discovered. But in our opinion, these results are unsuitable to Hamiltonian system due to its 
conservative property, here we only give a qualitative analysis for this phenomenon without 
proofs. When system is MS and within chaotic state, the coupling is stronger and trend 
to lock phases, this is exhibited by the platforms in Fig. 6(a), but due to the stochastic 
property of chaos, made the locking unstable, this is exhibited by the phase slips. Because 
the first possible locking is near vr, and sometimes tt locking is unstable, made ip stays a long 
time within (0, vr) (phase locking), then broke suddenly because the stochastic property of 
chaos, after 27r increase, (p travels again near vr, this time it has two choices, continue this 
slip if its still local unstable, or stay on locking state except local stable be changed. In our 



studies, we found that, in chaotic state, although phase locking doesn't exist, but still 
owns the trend to be vr locked as e increases, and the distribution for Ag also incline toward 
long distance, these properties made e decreased in chaotic state with e increases. 

In conclusion, we have investigated the way from non-MS to MS in coupled (j)^ lattices 
both for period and chaos. Due to Hamiltinian system's conservative property, exactly phase 
locking, AO = , is unreachable, but we found that, in periodic state, Lp trend to vr as e 
approach ^c, after e > Ec, A9 is oscillate around zero and the related amplitude decreases 
with e increases, on the critical point Ec, energy e own its minimum, also, different scalings 
has been found to describe these trends. MS is found to be connected closely with chaos in 
this model, random walk like phase locking is discovered, after analysed system's dynamical 
behaviors, we gave the qualitative explaining for this phenomenon. We argued that in some 
models, like 0^, Lyapunov exponent can be a suitable parameter for MS. But there are still 
some works ahead waiting for future research. The first work is to explain the phase locking 
and slip in chaotic MS strictly, we noticed that a recent work about phase synchronization 
also shown a random walk like phase locking |TE[ , its gave an example of nontrivial interaction 
between microscales and macroscales, so another work is to investigate the scalings for 
phase locking in chaotic states, also, to study these phenomena in large systems in also 
a interesting work. We wish this paper could deepen people's understanding of MS, and 
the observed phenomenon could potentially be investigated in many other related research 
in coupled Hamiltonian systems, e.g., elliptical galaxies, Hamiltonian version Kuramoto's 
model, transitions to equipartition in coupled large Hamiltonian systems, and we also hope 
this work could helpful in the research of quantum chaos. 

This research was supported by the National Natural Science Foundation of China, the 
Nonlinear Science Project of China, and the Foundation of Doctoral training of Educational 
Bureau of China. 
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Captions of figures 

2 2 

Fig. 1 Tfie time average energy for each lattice = ^im ^ !q + '^)d't-, i = 2 versus 
coupling strength s, critical coupling for MS is Sc = 0.0032. Total energy e = ei + 62 is 
plotted on the left vertical axis (denoted by hollow circle), on the critical point, e has its 
minimum, before and above £c, e decreases and increases monotonically. 

Fig. 2 With the initial condition Eqs. 4, when e = 0.002 < Sc, (a) plot the phase 
difference versus t, the slope rate 77 = 0.1, (b) plot phase orbits for two lattices, it is shown 
that two orbits moving in different domains without overlap, (c) and (d) corresponding to 
£ = 0.0035 > Sc, A9 oscillate with frequency / = 2.12, and two orbits be MS. (e) plote r] 
versus £ for £ < Sc and (f) plot / versus £ for £ > Sc- 

Fig. 3 (a) Probabihty distributions of cp for three different £S, circle for e — 0.001, upward 
triangle for s = 0.002 and downward triangle for e = 0.025, all below Sc- (b) Probability 
distributions for distance difference Aq, the couplings same to (a). 

Fig. 4 (a) The tt distribution of ip versus coupling for e < after logarithm, all data he 
on one line, the related scaling reads = ^^^^ with 7 = 1.96. (b) Oscillation amplitude 
A versus £, scaling relation reads A — A0 + Be~^ with a = 1 x 10~^. 

Fig. 5 Initial condition as Eqs. 7. (a) The summational energy e and each lattice's 
energy versus £, hollow circle represents the summation, (b) e and the largest Lyapunov 
exponent versus e. It can be found that in chaotic MS regions the summation decreases as 
£ increases. 

Fig. 6 £ = 0.025, the time evolution of phase locking in chaotic MS, A = 0.08. The inset 
window is enlargement of the rectangle region in original figure, phase locking and slips can 
be found. 
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